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Abstract 

We define a fixed point topological charge for the two-dimensional 
0(3) lattice (T-model which is free of topological defects. We use this 
operator in combination with the fixed point action to measure the topo- 
logical susceptibility for a wide range of correlation lengths. The results 
strongly suggest that it is not a physical quantity in this model. The 
procedure, however, can be applied to other asymptotically free theories 
as well. 
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1 Introduction 



Topological effects play an important role in the dynamics of asymptotically 
free field theories. In QCD instantons may be responsible for breaking the 
axial symmetry resolving the so-called U(l) problem In a large Nc limit 
the topological susceptibility relates the masses of the pseudo-scalars 77, rj' and 
K §. 

The topological susceptibility xt may be defined as the infinite volume limit 

x!'^^. (1) 

where Q is the topological charge and V is the space-time volume. In the 
two-dimensional 0(3) non-linear a-model it is a dimension two quantity that 
vanishes to all orders in the weak coupling expansion. From the perturbative 
renormalization group (RG) it is expected to scale according to the two loop (3 
function 

Xt a/3'exp(-4^/3), (/3^^). (2) 



It is a non-trivial task to recover the correct continuum results from lattice 
Monte Carlo simulations. A lattice topological charge definition is needed which 
returns even for large fluctuations reliable results. 

A 'geometric' definition proposed by Berg and Liischer 1^ is based on adding 
up the area of spherical triangles which are defined by the spin vectors in an 
elementary plaquette. As the contributions from all plaquettes are summed up, 
the internal space — the sphere described by the spin variables — is covered and 
if periodic boundary conditions are used one obtains an integer charge signify- 
ing the number of times this sphere is 'wrapped'. The topological susceptibility 
evaluated with this charge definition (and the standard action) completely failed 
to scale [p™. The reason was ascribed to special configurations called 'dislo- 
cations' M, which are dominant in the statistical average. Dislocations are 
non-zero charged configurations whose contributions to the topological charge 
come entirely from small localized regions where they become 'singular'. If 
the minimal action of dislocations is smaller than the continuum value of a 
one-instanton configuration (i.e. An) then dislocations will dominate the path 
integral and spoil the scaling behaviour, eq. (||) [|[||. 

Another definition goes back to DiVecchia et al. joj — for a recent discussion 
including the fixed point (FP) action see ref. Q. It is called field theoretical or 
plaquette definition and uses a 'naive' discretization of the continuum charge op- 
erator. This prescription does not yield integer values and to obtain continuum 
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results renormalization factors are needed. For large /3 these factors can be de- 
termined perturbatively, but for intermediate /3 one has to use non-perturbative 
techniques [p|-pl| . Results obtained with the field theoretical charge indicate for 
the susceptibility a behaviour consistent with scaling ^ ^ . 

A serious problem in these approaches is the role of the lattice artifacts, 
sensitive both to the form of the lattice action and the choice of the topological 
charge. A recent work suggests to use the FP action of a renormalization 
group transformation to study topological effects. In particular, an important 
feature is that the FP action has scale invariant instanton solutions (with an 
action value exactly An), and hence — as will be discussed in this paper — one 
can define a topological charge with no lattice defects. In ref. |]l3| and here the 
0(3) (T-model is considered, but the methods apply to other asymptotically free 
theories as well. The SU(3) gauge theory has been studied in refs. 14, 15|. A 



subsequent paper by one of us [ [17[ , will deal with the application of these ideas 
to CP""^ models and in particular to the CP"^ model. Some of our results were 



already presented in ref. |16 



The paper is organized as follows: First we review some results derived in 
ref. and define the FP field operator. This is followed by a closer look at 
instantons in the continuum, in a finite periodic volume and finally on a lattice 
using the FP action. We then define the FP topological charge and present 
some numerical results on classical solutions. In the last section we analyze the 
topological susceptibility evaluated in a Monte Carlo simulation. After a brief 
description of the methods used, we present the results which are followed by a 
conclusion and an outlook. 



2 RG results at the classical level 



2.1 Review of the RG transformation and its fixed point 



Let us briefly summarize some RG results which were developed in a previous 
paper jlSj. For a detailed discussion we refer the reader to this paper. 

We consider the 0(3) non-linear cr~model in two-dimensional Euclidean 
space defined on a square lattice. The partition function reads as follows 

Z= /oSe^'^-^t^). (3) 
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Here DS is the 0(3) invariant measure 

n 

and /3A{S) is a regularization of the continuum action 

pA^^U^) = f^ Jd^xd^Six)d^S{x), where S''{x) = l. (5) 

We perform exact RG transformations by a KadanofF type of bfocking, i.e. 
we divide the lattice into 2x2 bfocks labeled by indices tib- To each block we 
define a block spin variable R„b which is some mean of the spin variables S„ in 
the block. The block spins R„b form a lattice whose spacing is twice as large as 
the original one. An effective action is defined by integration over the original 
lattice: 

g-/3M'(R) ^ yj3Se-'3['4(S)+T(R,S)]^ (g) 

where T is the kernel of the RG transformation and its normalization 

^DRe-''^(^'^) = 1 (7) 



ensures the invariance of the partition function under this transformation. In 
the classical limit (/3 — > 00) the path integral is dominated by its saddle point: 

yl'(R) =min{^(S)+r(R,S)}. (8) 

The transformation kernel used in ref. has in the limit (3 ^ 00 the simple 
form: 



T(R,S)-«^( 



R«B • E • (9) 



nenB 



Here k is a free positive parameter of the RG transformation, which is tuned to 
make the FP action as compact as possible. As indicated by the free field theory 
in one dimension, the choice n = 2 gives the most short-ranged FP action. A 
fixed point of the transformation satisfies the equation 

ytpp(R) = min{^pp(S) + T(R, S)} . (10) 

This equation — called FP equation — fixes for arbitrary configurations {R} 
the value of the FP action. Starting from a lattice regularization of the con- 
tinuum action repeated RG transformations will drive the effective action to its 
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Figure 1: A multigrid is obtained by iterating the FP equation. 



fixed point. This takes on the form of a minimization in a multigrid of lattice 
configurations: 

^('=)(R)= min (yt(°)(SW) + r(SW,S('=-i)) + ... + T(R,S(i))].. 

{s<i),s(2),...,s('=)} L ' ' ' ' 'J 

(11) 

On each successive level — see figure ^ — the spin configurations become 
smoother, not only because the lattice spacing is halved, but also because the 
minimization tends to smooth out the fluctuations around a solution to the 
equations of motion. Hence one may choose for the action ^'"^ (S*^*^)) on the 
finest configuration {S*^*^^} any lattice discretization of the continuum action. 
The FP action is then obtained as the limit of A; cx) of ^'^'^■'(R). For 
practical purposes, however, only a few levels are needed and starting from the 
standard action on the lowest level the FP value is reached soon. 



2.2 Parametrization of the FP action 

In principle the above multigrid approach can be used to evaluate the FP ac- 
tion for arbitrary configurations to any precision desired. For practical calcu- 
lations, however, a parametrization of the FP action is needed. In ref. [p^ a 
parametrization has been obtained by fitting the known values of the action 
for ^ 500 configurations. That parametrization represented well the FP action 
on those configurations. However, to control the topological effects better, we 
decided to improve the parametrization further by including some small size 
topological solutions in the fitting procedure. We used several two-instanton 
solutions of the lattice FP action. While improving the fit for these instanton 
solutions, the new parametrization does not affect the quality of the fit for the 
previous configurations. 
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# 


type 


coupling 


type 


coupling 


type 


coupling 


1 


• — • 


0.61884 


•=• 


-0.04957 




-0.00932 


4 




0.19058 


/ 


-0.02212 


/ 


-0.00746 


7 


z 


0.01881 


z 


-0.00180 


X 


0.00658 


10 


J 


0.02155 


J 


0.00536 


J 


-0.00081 


13 


A 


0.00941 


A 


0.00488 


A 


-0.00225 


16 


X 


0.01209 


X 


0.00534 




0.00066 


19 




-0.00258 




-0.00173 




0.00146 


22 


E 


-0.01040 


E 


-0.00218 


E 


0.02720 



Table 1: Couplings used for the parametrization of the FP action including 
instanton configurations. The graphical notation is explained in the text. 



The resulting couplings are given in table ^ — together with a graphical 
notation of the corresponding operators. Let us explain here again the meaning 
of this notation. The parametrization of the action has the form 

^Fp(S) = ^ coupUng X products of ^"^n.^nj j (12) 

where is the angle between the two spins S„. and S„^- . Two dots connected 

with a line • • represent a factor in the action and the positions of the 

dots represent the lattice sites ni and rij respectively. Double, triple connected 
dots stand for the square, cube of the above factors. The operator, finally, is 
the product of all the factors 5??^. „^ as indicated by the lines in the figure. 
The quadratic and quartic couplings # 1,2,4,5,7,10,16 and 19 are determined 
analytically p3| , the others were determined with a numerical fitting procedure 
with the new instanton configurations added. 



2.3 The fixed point field 



As we shall see in section 3.2, FP operators are closely related to the FP field. 
The FP field is the fine field S'-'^-' in the multigrid solution of the iterated FP 
equation (0) as k goes to infinity. If the functional dependence of the solution 



on the first fine level S*^"'^^ on R is known, the FP field can be evaluated by 
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iteration. Below we construct the operator S^^^ = S(^'(R). (In the hmit k oo 
the solution {S^^^} of the iterated FP equation is identical to the solution {S} 
of the FP equation.) 



For smooth fields a quadratic approximation of the FP equation can be made 
which can be solved analytically. Consider a smooth configuration {R} where 
the spins fluctuate around the first axis: 



R„ 



v/1 



2 

riB 



Xn-B 



(13) 



where 2 components, and IxmbI ^ 1- Foi' the minimizing fine field we 

can make the ansatz 



c(i) 



(14) 



Inserted in the FP equation, the above expansion leads in leading order to the 
free field case. This was solved by Bell and Wilson ||l^ — for a brief review see 
for instance ref. — here we only report the relation between ??„ and Xn-B ■ 



^a(n, 



(15) 



Here a is given by 

a{n, Hb) = 



(27r)2 



-2iqj 



PFp(g) fj^ i-e 



(16) 



where p^piq) is the coefficient in the quadratic part of the FP action (given by 
the free field case). By iterating eq. (^5|) one can obtain the FP field operator 
in the free field case. The form is very similar to the above result, but with 
slightly modified parameters [Q. 



On smooth configurations {R} eqs. (|13H15|) give a good approximation. 
However we want to evaluate the fine field using a parametrization that per- 
forms well not only for smooth configurations {R}, but also for general ones. 
Our experience with the parametrization of the FP action suggests the ansatz 



S„=AA 



^a(n,nB)R„B + ^ /3(n, Ub, m^, 771^)^19^^ „j^R, 



(17) 



where A/" is a normalizing factor which ensures S'^ = 1. Like for the action it 
proved to be useful to use instead of the scalar product (1 — R,„^R„i^) between 
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the coarse spins at sites TOb and respectively, the angle ^, . In order 

to determine the coefficients /3, we numerically minimized the FP equation (|o|) 
using 60 configurations with lattice size 5 as input and stored the resulting fine 
lattices. The coefficients were then determined by minimizing the difference 
between the minimized fine spins and the parametrization dl] 



The numerical values of the coefficients a and (i are given in table g together 
with a symbolic notation of the corresponding operators. We chose a set of 
23 operators mainly because of their compactness. # 1-6 are the analytically 
determined coefficients a and # 7-23 are the numerically determined coefficients 
/3. The meaning of the graphical notation of the operators is the following: The 
dashed lines represent a 3 x 3 section of the coarse lattice grid. The cross + 



inbetween indicates the position n of the fine spin S„ in equation (17). The little 
square □ denotes the position n-g, of the coarse spin Rng . The two connected 

dots • • are the positions and of the spins whose angle \'Q'^^ ^, enters 

into the parametrization. Graphs obtained by trivial symmetry transformations 
are not drawn separately. 



3 Inst ant ons 



3.1 Instantons on the lattice 



In infinite volume continuum field theory topology is a well defined concept. 
Field configurations can be classified in topological sectors according to a 'wind- 
ing number' or topological charge In a lattice formulation this concept 
breaks down. When discretizing a theory, continuity in coordinate and internal 
space is lost. On the other hand topology is based on continuous transforma- 
tions of mappings which are separable in classes. In a discretized theory every 
field configuration can be continuously transformed into any other. If lattice 
configurations are sufficiently smooth, an unambiguous topological charge may 
be assigned. Conversely, for field configurations containing large fluctuations an 
interpolation is not unique and a charge definition becomes ambiguous. 

An additional problem arises due to the discretization of the continuum 
action. While the continuum action possesses scale invariant instanton solutions, 
this is generally not true for discretized actions. In particular, starting with non- 
zero charged configurations Liischer Q found that one can continuously lower 
the standard lattice action to zero by a local minimization in the spin variables. 

In reference ||l^ it was suggested to investigate the above problems with the 
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# type coeff 



# type coeff 



# type coeff 



1 \"^"'\ 0.59497 



2 



0.15621 



3 ; +; ; 0.08300 



0.00942 



-0.00171 



-0.00668 



7 ; +; -0.01228 



-0.02004 



9 IV; -0.03832 



10 



it- 



-0.05095 



11 ; +: 



0.01475 



12 



f::l: 



-0.00586 



13 tT;""""! -0.00303 



14 IVi i -0.00123 



15 T +i i -0.00118 



16 ; 



-0.01596 



17 i +i i -0.00090 
di-- 



18 



0.00140 



19 [X+i i -0.00128 



20 



i:t:j 



0.00447 



[p.—.-. 



21 i + 



0.00199 



22 



0.00304 



23 



-0.00015 



Table 2: Coefficients of the parametrization of tlie fine field. 
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aid of renormalization group methods and it was shown that the FP action has 
scale invariant instanton solutions. In this section we continue along these lines 
and construct instanton solutions of the FP action. These in turn can be used 
to study the performance of a proposed improved topological charge definition. 
But first let us review some basic facts about instantons in the continuum p!9| . 

In an infinite volume configurations with a finite action play a special role: 
At 'infinity' all spin variables point in the same direction and the space IR^ 
can be compactified by stereographic projection into a sphere 5*^. A finite 
action configuration is thus a mapping of a 'coordinate' sphere onto the internal 
sphere = 1. Such mappings can be classified by homotopy classes, with an 
integer number — the topological charge Q — characterizing the sectors. While 
configurations from the same topological sector can be continuously deformed 
into each other, this is not true for configurations with a different charge. The 
charge Q is the number of times the internal sphere is wrapped as the coordinate 
sphere is traversed. It may be defined as the integral 

Q ^ _L [d\ e,,, S • (a^S X 9,S). (18) 
Svr J 

and it is related to the action by the inequality 

^>47r|g|. (19) 

If for a given configuration the equality is satisfied, the configuration minimizes 
the action and is therefore a solution of the equations of motion. 

Let us now turn to the theory in a finite volume. By demanding periodic 
boundary conditions 

S(xi + Lm, X2 + Ln) = S(a;i, X2), where m, n & 'K, (20) 

we define the theory on a square torus of size L. In a finite volume every field 
configuration has a finite action and due to the periodic boundary conditions 
an integer topological charge Q associated with it. 

We can now explicitly construct pure instanton or pure anti-instanton con- 
figurations with an action ^ = 47r|(3| [|o). We use the plane coordinates defined 
by the stereographic projection to describe the solutions: 

2 Ui 

^ TTRF 

(21) 

_ l-\u\^ 
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where u — ui + iu2 ■ The instanton solutions at the boundary of equation ( |19[ ) 
satisfy the Cauchy-Riemann equations for u being an analytic function in z = 

Xl + iX2- 

{di + id2)u = 0. (22) 

The solutions are doubly-periodic meromorphic functions called elliptic func- 
tions pO| . They can be written as 



- = ^n^fe5fy' -th ta. = j:k. (23) 

1=1 ^ i=l i=l 



Here the integer k is the topological charge of the solution and c, ai, . . . , Ofc 
and bi, . . . ,bk are complex numbers. (t(z) is the Weierstrass cr-function with 
half-periods lu — L/2 and u)' — iL/2. Using Cauchy's theorem one can show 
that there are no solutions with topological charge equal to one Hence we 
are forced to construct charge two instanton solutions. Specifically, we set k = 2 
and c = 1. A reasonable definition for the instanton size is 

p = imin{|ai - |ai - 62|}- (24) 



Let us now turn to the construction of instanton solutions on the lattice. 
As was pointed out in ref. ||l^ the FP action allows scale invariant instanton 
solutions. Since we use this fact to construct instantons on the lattice, it is 
appropriate to repeat the statement: 

If a given configuration {R} satisfies the equations of motion for the FP 
action and it is a minimum, then the solution {S} of the FP equation ( p!o| ) 
satisfies the equations of motion also. Moreover, both configurations yield the 
same value for the action. 



The proof is quite simple: If {R} is at a local minimum of the FP action 
.App, variations with respect to {R} will vanish: 



JAp(R) 

S'R-n-B 



1^ I ^ ^ S„ + Rng I -f^iiB ' ^ ^ ' 



S„ = 0. (25) 



Here S = S(R) is the solution of the FP equation with {R} as coarse input 
configuration. Since {R} is a local minimum eq. ( p5|) implies 

R„3^ (26) 

Consequently, we have 

T(R, S) 0. (27) 
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Since the transformation kernel T(R, S) > 0, the configuration {S} gives its 
minimum for fixed {R}. Because the configuration {S} minimizes the right 
hand side of the FP equation, it minimizes the FP action y^Fp(S) separately. 
Therefore it is a solution of the equations of motion. Since T(R, S) = 0, both 
actions have the same value: ylFp(R) — ^fp(S). This concludes the proof. 

The reverse of the above statement is in general not true for arbitrary config- 
urations. Although a fine configuration which is a solution of the FP equations 
of motion and which is used to construct a coarse configuration by means of 
the above blocking equation (p6|), locally minimizes the right hand side of the 
FP equation, this minimum does not necessarily coincide with the absolute one. 
This in fact prevents the existence of arbitrary small instanton solutions. 

Using the above ideas it is clear how to construct instanton solutions of 
vApp. We naively discretize the continuum instanton solution on a very fine 
lattice with spacing qq = 2~^a. After performing k blocking steps as defined 
by equation (p6|), we obtain a configuration {R} on a lattice with spacing a. 
In the limit k oo we recover the continuum solution, which is of course 
a local minimum of the continuum action. Since all the successive blocked 
configurations minimize the transformation kernels in the iterated FP equation, 
the configuration {R} is a good candidate for a lattice solution of the FP action]^. 
We may solve the FP equation for {R} to check whether it is still a solution. If 
{R} is a solution, then the multigrid minimization procedure should lead to the 
same configurations on finer lattices as those which were used in constructing 
{R} by blocking. 



3.2 Definition of the topological charge on the lattice 

In the following we define a topological charge operator based on the multigrid 
solution of the FP action. We evaluate the FP topological charge by means of 
the solutions of the FP equation (p^. Under a RG transformation an operator 
0(8) transforms into C(R) on the coarse lattice as 

0'(R)e~'5'-^'W = o(s)e-'3[-4(s)+T(R.s)]^ (38) 

In the limit /? — > cxd the path integral on the right hand side is approximated 
by its saddle point and we obtain 

C(R) = 0(S(R)), (29) 

^ It will be a solution, unless the size of the instantons is too small with respect to the 
lattice spacing. 
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where the spin configuration {S(R)} is the solution of the FP equation (p^. 
Repeated application of this transformation will single out the operator with 
largest eigenvalue. Since the topological charge is expected to be a marginal 
operator, we may obtain it as the limit 

Qpp(R) = lim Q(S('=)(R)). (30) 

k — ^oo 

Here Q is some standard lattice charge definition and {S*^*^)} is the solution of 
the iterated FP equation ( pl| ) on the lowest level in a fc level multigrid (see 
figure |l|). In other words, the FP topological charge is a standard topological 
charge evaluated on the FP field. 

Note that S*^''^ (R) becomes increasingly smooth as k grows: first, because the 
corresponding lattice spacing ao — 2^^a decreases, secondly, because S*^'^^(R) 
becomes almost a solution to the equations of motion. Consequently, any sensi- 
ble definition of the topological charge can be used on these configurations, the 
final result will not depend on this choice. Nevertheless, it is more convenient 
to use the geometric definition since it is stable against small variations of the 
field and hence k, the number of levels in the multigrid minimization could be 
kept small with this definition. For a review of the geometric definition we refer 
the reader to ref. 

One can easily show that with this definition of the topological charge there 
are no dangerous dislocations. More precisely, one has 

^PP(R) >47r|Qpp(R)|, (31) 

for arbitrary configuration {R}. The corresponding statement is true in the 
continuum, hence it is also true for S'^'^''(R) for k ^ oo. Eq. (^ij) follows then 
by observing that the contribution of the T terms in eq. ( |lO| ) is non-negative. 
We are now ready to discuss the numerical aspects of classical solutions. 



3.3 Classical numerical results 



Following the above program we naively discretize two-instanton solutions on 
the torus of various sizes on very fine lattices. We found that four to five block- 
ing steps are sufficient to make any lattice artifacts of the original discretization 
negligible. On the finally blocked configurations we can measure several quan- 
tities. On the coarse configuration itself we measure the standard action and 



the parametrization of the FP action presented in section 2.2. Performing a 



minimization on a multigrid with three finer levels we measure the exact FP 
action and on the finest level the FP charge. Using the instanton radius given 



13 



by eq. (^4|), we get a parameter which characterizes weh the breakdown of the 
blocking to obtain instanton solutions on coarse lattices: In figure |^ it can be 
clearly seen that below an instanton radius of p < 0.7a the configurations are no 
longer instanton solutions and we shall say that the instanton falls through the 
lattice. It is gratifying to see that the FP charge immediately falls off to zero, as 
the FP action drops below the continuum value. Furthermore, figure |^ demon- 
strates how well the parametrization for the FP action is suited for instanton 
configurations. The deviation from the exact value is quite small, in particular 
the parametrized FP action is only marginally smaller than the continuum value 
in the region above the point where the instanton falls through the lattice. In 
contrast, the values of the standard action are quite different from the contin- 
uum ones. These numerical findings support the above statement, that there 
are no dangerous dislocations present when using the FP action together with 
the FP topological charge. 



2.2 




0,6 - A standard action 

• parametrized FP action 

X minimized FP action 
0.2 FP topological charge 

0.0 ' 

0.4 0.5 0.6 0.7 0.8 0.9 

p/a 

Figure 2: Actions and charge of instantons with radii of the order of one lattice 
spacing 



In numerical simulations we use a parametrized form of the FP action and 
the parametrization in eq. ( |l7| ) for the FP field instead of the time consum- 
ing minimization procedure. One then is interested if the use of these two 
parametrizations has an influence on the existence of dislocations. The curves 
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referring to the FP action and FP charge in fig. |^ have a non-analytic break 
at p/a ~ 0.7. Our parametrization does not fully reproduce this behaviour, 
and so one might expect that observable deviations will occur in this instanton 
region. We systematically searched for minimal action configurations with a 
parametrized action lower than the continuum value in the Qfp'' — 2 sector. 
We found a minimal action of 1.84 • 47r (compare this with the value 0.93 • Air we 
found using the standard action and geometric charge). This value can be as- 
cribed to the not exactly accurate parametrization of the FP field. If we actually 
solve the FP equation ( |l0[ ) for this 'dislocation configuration' we get the correct 
charge Q — Q. Nevertheless, as one has to use parametrizations for Monte Carlo 
simulations, such configurations could be dangerous. On the other hand, the 
search for dislocations reveals the weakest point of the parametrization which 
performs very well in other cases (c.f. figure ||). What actually counts is not how 
the parametrization works for some configurations that were specially sought for 
their bad performance, but how well it performs for configurations in thermal 
equilibrium occurring in a Monte Carlo simulation. The results of a test of 



this performance — presented in section |4T| — shows that indeed there is no 
problem. 

We also searched for the minimal action configuration in the Q = \ sector 
and did not find a configuration with an action below the continuum action. 
This is not astonishing as there are no one-instanton solutions on a torus. 



4 Topological Susceptibility 



If the topological susceptibility is a well defined physical quantity that is renor- 
malization group invariant, then one expects that it scales like a (mass)^ in the 
continuum limit. One additionally measures a second quantity, e.g. the corre- 
lation length ^ and builds the dimensionless product xt S,^ which should go to 
a constant in the limit ^ ^ oo. Earlier Monte Carlo calculations do not show 
convincingly whether this is the case. Furthermore, perturbative considerations 
indicate that in the 0(3) model there might be a problem with the topological 
susceptibility in the continuum limit. 

One may calculate the contribution of instantons in the continuum using 
a semiclassical expansion. The probability density to find an instanton with 
topological charge Q ~ \ and size p <C 1/A is pl|,p^: 

Pi - 'A (32) 
P 

where A is the scale parameter of the model. Using the renormalization group 
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one can show, that eq. (^2|) is exact as p —^ assuming the smaU instantons 
form a dilute gas. The infrared divergence at p = in eq. ( p2[ ) indicates, how- 
ever, that this assumption is not true: small instantons are not suppressed, but 
contribute strongly to the susceptibility. Such a dominance of small instantons 
is also indicated by numerical studies, trying to determine the instanton size 
distribution with different methods |pl] , [l^ . 

On the lattice, however, not the whole range of the instanton radius is 
probed. The lattice cuts the contribution of small instantons because there 
is a smallest possible radius before the instanton falls through the lattice (c.f. 
figure H). Hence the measured topological susceptibility is always finite, ft is not 
excluded, however, that it raises boundlessly with increasing correlation length 
as the perturbative considerations suggest. 



4.1 Numerical Results 



Using the perfect lattice action and the perfect charge we performed extensive 
Monte Carlo simulations at correlation lengths in the range E {2 — 60). In 
order to avoid finite size effects we kept the ratio L/S^k, Q constant. The corre- 
lation length was obtained from the long distance fall off of the zero momentum 
correlation function. 

We determined the topological charge using both the geometric definition 



and the definition of the FP charge given in section 3.2. For the measurement 
of the FP charge we used the geometric definition of the charge on a finer lattice 
of the multigrid with the Monte Carlo generated lattice as coarsest level. To 
determine the configuration S(R) on the fine lattice one can either minimize 
the FP equation (which is very time consuming) or use the parametrization of 
the dependence on R given in section |2.3| . We denote the corresponding charges 
Qcoarsc for thc gcomctric charge, Qpar. i. lovoi for the charge measured on the first 
finer level using the parametrization of the fine field etc. 

For two /? values we compared the results of using the parametrization on a 
finer level and of minimizing on a multigrid. As is shown in table ^, the results 
were found to be consistent within the statistical errors. This confirms, that the 
parametrization performs well for configurations occurring in a Monte Carlo 



simulation (c.f. the discussion at the end of section |3.3| ). In order to test if the 
result of going to a lower level is already stable, we also calculated the charge 
for a few (3 values on the second finer lattice. As also can be seen in table ||, 
the values on the first and the second finer level were found to be consistent 
within thc statistical errors. Thus it is sufficient to calculate the fine field only 
on the first finer level. This is not astonishing as the maximal angle between two 
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/3 


(Qcoarsc) 


(^par. 1. level) 


(Qpar. 2. level) 


{Qra'm. 1. level) 


(*3miii. 2. level) 


0.6 


2.56(5) 


1.96(4) 


1.93(4) 






0.7 


3.39(4) 


2.48(3) 


2.45(3) 


2.45(3) 


2.45(3) 


0.85 


5.62(11) 


4.31(9) 


4.30(9) 






1.0 


6.38(8) 


4.86(6) 


4.82(6) 






1.0 


6.33(6) 


4.78(5) 




4.75(5) 





Table 3: Results of test MC Simulations in the 0(3) model indicating that it 
is sufficient to measure the charge only on the first finer parametrized level. 
Minimizing or going to a lower level yields the same result. 



/3 


L 




(^coarse) 


Xcoarse ' ^' 


(*3l. level) 


Xl. level ■ 


0.51 


10 


1.6049(42) 


1.87(1) 


0.0482(4) 


1.329(8) 


0.0342(3) 


0.6 


14 


2.1960(46) 


2.57(2) 


0.0631(6) 


1.90(2) 


0.0467(4) 


0.685 


20 


3.012(14) 


3.57(3) 


0.0810(10) 


2.65(2) 


0.0601(7) 


0.7 


20 


3.186(15) 


3.36(3) 


0.0852(11) 


2.52(2) 


0.0640(8) 


0.85 


40 


6.057(17) 


5.73(4) 


0.1314(12) 


4.38(3) 


0.1004(9) 


1.0 


70 


12.156(34) 


6.36(5) 


0.1918(19) 


4.80(4) 


0.1448(15) 


1.1 


130 


20.397(86) 


10.04(12) 


0.2472(36) 


7.69(9) 


0.1893(27) 


1.2 


180 


34.44(30) 


8.23(14) 


0.3013(73) 


6.11(10) 


0.2237(53) 


1.3 


340 


58.06(37) 


12.01(29) 


0.3502(96) 


8.94(22) 


0.2607(72) 



Table 4: Results of MC Simulations with the FP action. 



neighbouring spins halves as one goes one step down to a finer level. So even on 
the first finer level the maximal possible angle between two neighbouring spins 
is 90° and there is practically no ambiguity left for the topological charge. 

In table ^ we report the results of the simulations for the correlation length 
^, the geometric charge (Qcoarse) the FP charge {Ql i^^^^) evaluated on the 
first parametrized finer level. Using the geometric charge as well as the FP 
charge, we build the dimensionless quantity x* • of topological susceptibility 
and correlation length. 

Figure ^ shows the results for the topological susceptibility. Clearly, no 
scaling is seen even at correlation lengths as large as 60. Both curves, the one for 
the geometric charge and the one for the FP charge, are rising and no flattening 
occurs at the largest correlation lengths. There is a significant difference between 
the topological susceptibility built with the geometric charge and the one with 
the FP charge. The value of the geometric charge lies several standard deviations 
above the value of the FP charge. Furthermore, the difference is slowly growing 
with increasing correlation length. 



17 



0.40 



• coarse level 

0.30 - ♦ finer levels I 



* 0.20 



+ 

0.10 - + 
0.00 I ' ' ' — ^ 



+ 



I 



I 



10 100 

Figure 3: Results of Monte Carlo measurements of the topological susceptibility 
at different correlation lengths. 



5 Conclusion and Outlook 



The definition of the FP topological charge is based on the FP field operator. 
The FP field can be evaluated to any precision desired by solving the iterated 
FP equation on a k level multigrid. As was demonstrated by fig. ^ and can be 
shown analytically for sufficient large multigrids the FP charge together with 
the FP action has no lattice defects whatsoever. 

For use in MC simulations however a parametrization of the action and 
the charge are needed. Instead of parametrizing the charge directly, we have 
parametrized the solution of the FP equation which is to be iterated to obtain 
the FP field. The accuracy of the parametrizations have been rechecked in MC 
simulations. 

The partition function of the lattice tr-model is — as indicated by a semi- 
classical approximation [pT| — dominated by small sized topological excitations. 
These unphysical fluctuations are the cause for the divergence which is seen in 
fig-l 
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Although we studied the cr-model in this paper, the methods derived are 
quite general and can be applied to other asymptotically free theories such as 
CP"^"^ models or SU(N) gauge models. For instance, in CP" models a lowest 
order semiclassical approximation estimates for instanton contributions pl|p2| 

/~y" dpp""-^. (33) 

If N > 4, the contribution of short distance topological excitations is small rela- 
tive to that of the physical ones and we expect to see a scaling of the topological 
susceptibility according to the perturbative RG. In fact in the CP^ model, which 
is studied by one of us ||l^ , the topological susceptibility already exhibits the 
expected scaling behaviour. 
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